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ABSTRACT 

The planet formation process and subsequent planet migration may lead to configurations resulting 
in strong dynamical interactions among the various planets. Well-studied possible outcomes include 
collisions between planets, scattering events that eject one or more of the planets, and a collision of 
one or more of the planets with the parent star. In this work we consider one other possibility that 
has seemingly been overlooked in the various scattering calculations presented in the literature: the 
tidal capture of two planets which leads to the formation of a binary planet (or binary brown dwarf) 
in orbit about the parent star. We carry out extensive numerical simulations of such dynamical and 
tidal interactions to explore the parameter space for the formation of such binary planets. We show 
that tidal formation of binary planets is possible for typical planet masses and distances from the 
host star. The detection (or lack thereof) of planet-planet binaries can thus be used to constrain the 
properties of planetary systems, including their mutual spacing during formation, and the fraction of 
close planets in very eccentric orbits which are believed to form by a closely related process. 
Subject headings: accretion, accretion disks — planets and satellites: general — celestial mechanics 
— methods: TV-body simulations — planets and satellites: formation — stars: 
binaries: eclipsing — stars: low-mass, brown dwarfs — stars: planetary systems: 
formation — stars: planetary systems: protoplanetary disks 



1. INTRODUCTION 

The Earth-Moon and the Charon-Pluto systems are 
sometimes referred to as double (or binary) planets, i.e., 
binary systems consisting of two planets whose center 
of mass orbits a central star. These rocky systems are 
likely to have formed by the fi ssioning of a more massive 



planet due to a giant impact ( Ha rtmann fc Davis 1975 
Lin|1981 [Canup fc Asphaug|20IIT ) . In this study, we are 
interested in binaries ot gas giant p 



planets (or even brown 
dwarfs) for which a fission origin is unlikely. Specifi- 
cally, we investigate whether such systems can form by 
tidal interactions and their implications for planet forma- 
tion. The discovery of such systems in current or future 
searches fo r planet transits (suc h as the Kepler |Basri 



et al.|[2005 and CoRoT missions [Auvergne et al.||2UU9 



IS exciting for a number of reasons, i'hese include (i 
the fact that the existence of binary giant planets could 
provide strong observational evidence for tidal capture 
as a viable astrophysical mechanism, (ii) Binary plan- 
ets would allow for new and important tests and models 
of planetary dynamics early in the formation process of 
planetary systems, including their mutual spacing dur- 
ing formation, (iii) They would also allow for studies 
of long-term planetary dynamics, including current mea- 
sures of internal structure via apsidal motion, and spin- 
orbit interactions, (iv) Binary planets would provide un- 
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precedented accuracy for determining masses, radii, in- 
ternal structure, etc. (v) It seems quite possible that an 
eclipsing set of planets, transiting a parent star, would 
provide more information on the oblateness of the plan- 
ets than a s imple transit hght curve (see, e.g.. Carter fc 



Winn 2010), especially given that such a system would 
be ex]3ected to be rotating more rapidly than a single 
planet orbiting the parent star, (vi) Finally, if there is 
any significant amount of "magnetic braking" in a close 
binary planet system, the two planets could actually be 
driven into Roche lobe contact, leading to mass trans- 
fer between the planets. Such a planetary mass-transfer 
system would most likely be quite stable and very long 
lasting. 

The formation of a binary planet is intimately linked to 
the evolution of the protoplanetary disk from which the 
planets have formed and within which they evolve. When 
the two planets are still embedded in a gaseous disk, they 
may migrate inwards or outwards, transferring orbital 
angular momentum to or from the disk. If the migration 
is relatively slow, the two planets may evolve into an 
isolated low-order mean-motion resonance and will then 



migrate together locked in this resonan ce (Lee fc Peale 



2002| Papaloizou & Szuszkiewicz 2005). Such a stable 
configuration precludes the formation of binary. How- 
ever, if sufficient eccentricity is somehow induced during 
the migration process, or if the migration process is suf- 
ficiently fast to push the system through the resonance, 
neighboring resonances can destabilize the system. More 
generally, the orbits of two planets become dynamically 
unstable when the fractional differe nce of their orbital 
radii becomes sufficiently small (Gladman 1993). The 
presence of a disk can inhibit the development ot an insta- 
bility by limiting the growth of eccentrities. But once the 
stabilizing influence of the disk disappears, e.g., because 
of a decrease of the disk mass (more specifically, the disk 
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surface density) or because the planets have grown suffi- 
ciently by accretion from the disk, the planet orbits can 
become unstable. This will generally lead to a dynam- 
ical scattering event, in which the configuration o f the 
planets can cha nge drastically. Ford & Rasio f2008) and 
Chatterjee et al._ (2008) have systematically studied such 
scattering events m systems with two and three planets 
and found that the results of such scatterings could be: 
(a) the collision and merger of the two planets, (b) the 
collision and merger of a planet and the host star, (c) the 
ejection of one of the planets, or (d) a quasi-stable con- 
figuration in which both planets remain orbiting the host 
star after ~ 10^ orbits. However, these authors did not 
consider another possibility: the formation of a binary 
planet. The latter may result either from a thr ee-body 



exchange o r more simply from a tidal capture (Fabian 



et al. 



1975) 



A tidal capture occurs when the two planets get suffi- 
ciently close - typically within a few planet radii - but do 
not collide directly. Such close encounters induce tidal 
oscillations in one or both planets, converting orbital 
energy into oscillation energy (which is eventually dis- 
sipated as thermal energy) and leaving a bound binary 
planet (at least temporarily) in orbit about the host star. 
During subsequent periastron passages, energy may be 
exchanged either to or from the tides (unless the tidal 
oscillations have been completely damped in the mean 



time) , and the evolution is formally chaotic (Mardling 
1995 ) in much the same way as an unstable three- body 



orbit (again with the possibility of dissociation). Dur- 
ing this phase, the tides can be extremely large (because 
the tidal amplitude is additive) , and it is likely that non- 
linear fluid processes such as shocks operate to make en- 
ergy dissipation quite efficient, limiting the amplitude of 
the tides. However, once the system has dissipated a 
sufficient amount of tidal energy, it will cease to behave 
chaotically and will circularize in the normal fa shion with 



the t idal amplitude never being able to grow (Mardling 
1995 1. In the case that tidal capture occurs in the tidal 
field of a third body (in the present case, the star), the 
tidal energy exchange process can be doubly stochastic. 
However, to simplify the modeling process, we will as- 
sume in this work that all energy deposited in the tides 
is dissipated before the next periastron passage, consis- 
tent with severe tidal damping during this phase as dis- 
cussed above. The validity of this assumption depends 
on the efficiency of the tidal dissipation p rocess , often 



parametrize d using the tidal Q value (e.g., [Goldreich fc 
SotCTlp66| 



in a three-body exchange, a binary planet forms by 
the "exchange" of the outer planet into the "inner bi- 
nary" (the central star -|- inner planet) forming a new 
inner binary (the two planets). If there is no dissipa- 
tion of orbital energy, the newly formed pair is suscep- 
tible to dissociation because, unlike a normal two-body 
tidal capture, it is continuously forced by the tidal field 
of the star. The binary orbit will exchange energy and 
angular momentum with its center-of-mass orbit around 
the star in a random-walk fashion (the system is for- 



mally chaotic; Mardling|2008 1 , until sufficient energy has 
been transferred to the binary to dissociate it. This can 
only be avoided once enough orbital energy is perma- 
nently removed from the binary to make the system sta- 
ble against chaotic interactions. This can occur through 



(point- mass) interactions with other bodies (e.g., a third 
planet or a background of planetesimals), similar to the 
mechani sm proposed for the formation of Kuiper-Belt 
binaries ( Goldreich et al.|2002 ), or through tidal dissipa- 
tion if the planets are close enough. 

In both scenarios, tidal dissipation may play a key role 
in forming a binary planet, but, in the first case, the 
tidal capture has to operate in one (or possibly a few) 
encounters, while, in the second, tidal dissipation may 
operate over the much longer timescale of the transient 
binary state. We therefore expect that, in the first case, 
the post-capture orbital separation of the binary planet 
is a few planet radii, while, in the second, it could be as 
wide as is stable on a long time timescale. 

In this study, we consider primarily the mutual tidal 
capture of two planets in a "dynamically active" plane- 
tary system. We show that, for reasonable assumptions 
about the tidal coupling and dissipation of gas giant plan- 
ets, mutual tidal capture of planets into a planet-planet 
binary is a relatively generic feature of dynamically sys- 
tems. In this regard, it has recently been shown that the 
long-standing issue of the relatively high eccentricities 
of extrasolar planets can be understood if the configura- 
tion of the newly-formed planetary system after the gas 
disk has dissipated is dynamically active, so that planet- 
planet scattering can operate and eject planets whi le in- 
creasing the ec centricities of the r emainin g planets ( Ford 
& Rasio 2008; Juric fc Tremaine 2008; Chatterjee et a, 



2008|. It is not obvious that planetary systems forming 
via tne usual core accretion scenario should result in dy- 
namically active systems after the gas disk phase. How- 
ever, simulations starting with protoplanets in a gas disk 
and following both gas physics and A'^-body dynamics 
have show n that such res ulting confi gurations are indeed 
possible ( Thommes et al. 2008 , Jvlatsumura et al.pOlO) . 

Based on our calculations, we present the relative prob- 
abilities of the possible outcomes of dynamically active 
systems (including, e.g., planet-planet collisions, tidal 
capture, ejection). We discuss the potential observabil- 
ity of a planet binary in relation to the likelihood that 
planetary systems are dynamically active early in their 
lifetimes. 

In our study we also allow for the possibility of tidal 
capture of brown dwarfs (BDs) in close orbit around a 
hydrogen-burning star. Brown dwarfs probably form in 
a very different way from, and at a different distance 
than, gas giant planets. While gas giants likely form 
via core accretion, any brown dwarf forming out of a 
circumstellar disk will do so via gravitational fragmen- 
tation, and at fairly large distances from the host star 
(> 10 ^ AU) , where the disk is Toomre unstable ( Toomre 
1964[ ) and can cool su fficiently so that unstable clumps ^ 



can coll apse (see, e.g., Stamatellos fc Whitworth 2009 



also see Boss||1997|). We do not address the issue of how 



to get the brown dwarfs to distances of < 1 AU from the 
host star, but note that such systems are known to exist 
(e.g., CoRoT 15b contains a bro wn dwarf orbiting a n F 
star with an orbital period of 3 d; Bouchy et al.|2010 ). In 
any case, as we will show, if multiple brown dwarfs exist 
on such close orbits, they can easily become bound in a 
BD-BD binary via tidal capturejj 

^' Throughout this paper, we do not sharply distinguish between 
planets and brown dwarfs and often refer to all types of sub-stellar 
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In §2 we present the results of numerical simulations of 
dynamically active planet systems for a range of masses 
and distances from the parent star. In §3 we describe 
the detectability of such binary planets in radial velocity 
searches and in transit studies. 

2. NUMERICAL SIMULATIONS 
2.1. Numerical Method 

Since the dynamical processes leading to planet-planet 
tidal capture are the same as those that lead to planet- 
planet collisions, any numerical method which can prop- 
erly treat "dynamically active" planetary systems should 
serve as an appropriate base upon which we can add a 
treatment of tidal dissipation. We use the Fewbody inte- 
grator, which is designed for strong small- A'^-body gravi- 
tational encounters (Fregea u et al.|2004| ). We have tested 
that the code is suitable for simulating dynamically ac- 
tive planetary systems by compari ng our results to the 
simulations of Ford & Rasio ( 2008 1 (see Appendix) . 

In order to include tidal capture, we also had to in- 
clude a treatment of tidal dissipation in Fewbody. We 
use the functional fits to the energy dissipated by the 
close passage of two polytr opes as presented in [Porte- 1 
gies Zwart & Meinen ( 1993J ). When the two planets, of 
mass rrii and radius Ri, encounter each other with close 
approach distance less than 10(i?i -I- R2), we calculate 
the energy dissipated in the encounter by treating each 
planet as an n = 3/2 polytrope and reduce the planet- 
planet relative velocity (at pericenter) accordingly. We 
only include tidal dissipation in impulsive interactions, 
requiring - somewhat arbitrarily - that the close ap- 
proach distance is less than 10% of the previous local 
maximum in the planet-planet distanceF] This provides 
a lower limit on the overall tidal dissipation. If the close 
approach distance is less than i?i + R2 we assume the 
planets merge instantaneously, conserving mass and mo- 
mentum. 



2.2. Stability Criterion 

Fewbody automatically terminates a calculation when 
an unambiguous, dynamically stable configuration has 
been obtained. To accurately test for the dynamical sta- 
bility of a binary planet orbiting a star, w e use the newly 
improved algorithm of Mardling (2008), which utilizes 
the concept of orbital resonance overlap to test for chaos 
and, ultimately, instability. The resonance overlap algo- 
rithm in its current form is only approximate for small 
values of the outer eccentricity Cout, i-e., the eccentricity 
of the outermost body (in this context the parent star) 
about the center of mass of the innermost binary (the 
planet-planet binary). We thus set Cout ==0.1 in the al- 
gorithm when the outer eccentricity is smaller than this 
value. In addition, we also require that the apocenter of 
the planet-planet orbit is within 2/3 of its Hill radius at 
pericenter. Specifically, we require 
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Fig. 1. — Histogram of semi-major axes of circularized binary 
brown dwarfs resulting from tidal capture events in model t21 (with 
brown dwarfs masses of 30Afj and 70 Afj, respectively). 



where ajn and Cin are the semi-major axis and eccentric- 
ity of the newly formed planet-planet binary, and the ex- 
pression on the right is 2/3 the radius of the Hill sphere 
of the planet-planet pair evaluated at apocenter. 

2.3. Starting Conditions 

We adopt initial co nditions for our sim ulations that 
are similar to those of Ford & Rasio ( 2008 1 . The planet 
masses are either mi = 10"'^ Mq and m2 = 0.5 x 
10"^ Mq for the "gas giant" case, or mi = 70 x 10"^ Mq 
and 7712 = 30 X 10~^ Mq for the "brown dwarf" case. In 
all cases, we fix the central star's mass at Af* = IMq, 
and planet 1 initially orbits with semi-major axis ai.init < 
a2,init- We use three initial values for ai: ai.init — 0.2, 
1, and 5 AU. The initial semi-major axis of planet 2 is 
set so that the system is not Hill stable, with a2,init ran- 
domly distributed uniformly between 0.9ai^init(l -I- Ac) 
and ai.init (1 + A c), where A^ = 2.4(mi/M*-Fm2/M^)i/3 



objects collectively as planets. 

''' This means that no tidal damping is applied for eccentricity 
< 0.8 in a bound planet-planet system. 



(Gladman 1993). The planet initial eccentricities eiinit 
are randomly distributed uniformly between and 0.05, 
and the relative inclination of the planet orbits is dis- 
tributed uniformly between 0° and 2°. All remaining 
orbital elements are sampled uniformly in their allowed 
range. For reference, all model parameters are shown in 
Table m 

2.4. Thermal Bloating of the Interacting Planets 

In this study, we assume that each of the pair of plan- 
ets/brown dwarfs was born early in the evolutionary his- 
tory of the parent star. In particular, we adopt the hy- 
pothesis that the planets/brown dwarfs were either born 
separately in the unstable collapse of a massive accre- 
tion disk, or otherwise were driven to migration toward 
the parent star via an accretion disk. In either case, the 
planets or brown dwarfs were necessarily young when 
they captured each other, i.e., had an age of ~ 1 — 30 
Myr. Planets, and especially brown dwarfs, are quite 
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Fig. 2. — Histogram of semi-major axes of circularized binary 
planets resulting from tidal capture events in model tl6 (with 
planet masses of 1 Mj and 0.5 Mj, respectively). 

th ermally bloated at t hese young ages (see, e.g., Fig. 4 



of Nelson et al 
empirical htting 



of mass > 
follows: 



15 M, 



loos'). We have constructed a scmi- 
ormula for the radius of brown-dwarfs 
as a function of evolution time, i, as 



i?(t)=0.79VM,/A/jt//'i?j 



(2) 



where tg is the evolution time in units of 10^ yr. Mi is 
the mass of the planet, and Mj and Rj are the mass and 
radius of Jupiter, respectively. This expression works 
quite well for masses above ~ 20 Mj and for ages in the 
range of 0.1-30Myr. For any object whose radius falls 
below 1 i?j based on this expression, we simply fix the 
radius at 1 i?j |j 

Tidal dissipation is a strong function of the radius of 
the planet or brown dwarf compared to its separation 
from the object with which it is interacting. Hence, the 
inclusion of thermal bloating is potentially important, as 
it allows tidal capture to occur at larger initial separa- 
tions. 

After a planet-planet binary is formed via tidal capture 
and the resulting star-pl anet-plane t syste m is deemed 
dynamically stable by the Mardling (2008) stability cri- 
terion, we stop the calculation and record the semi-major 
axis, Oin, and eccentricity, ejn, of the planet-planet bi- 
nary. In post-processing we assume the orbit is quickly 
tidally circularized, and set Odrc = ain(l ^ e^J. Fig. UA 
shows a histogram of the circularized planet-planet semi- 
major axes resulting from tidal capture events in model 
t21 (for a pair of brown dwarfs). There is a clear peak 
just above Och^j ~ 2Rq, with a tail that extends out to 
~ 8Rq. Fig. l2]shows the distribution for model tl6 (for 
a pair of giant gas planets). 

2.5. Illustrative Scattering Results 



* Note that this prescription may underestimate the radii of 
very clos e Ju pitcr-mass planets that could be inflated due to tidal 
heating fBodonhcimer ot al. 2001). 
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Fig. 3. — Evolution of a planet/brown dwarf binary system from 
model t21 in the x—y plane at the start of the calculation (upper 
panel), and at the end of the calculation (lower panel). The star is 
represented by a black trail, and the planets by red and blue trails. 
The trails fade with time so that the darkest points along a curve 
are the most recent. Dissipative tidal encounters are shown as open 
circles at the point where they occur. The system becomes active 
quickly after the start of the calculation. Over time the planets 
exchange position relative to the star. The planets suffer a weak 
dissipative tidal encounter at late time but eventually collide and 
merge. 



To get a better feel for how the dynamics unfolds in 
a simulation resulting in tidal capture, we have plotted 
the evolution of the star and 'planet' positions for a set 
of representative simulations from model t2lP| Fig. [s] 
shows a typical simulation ending in a merger of the two 
planets. The system becomes active quickly after the 

^ Model t21 is a simulation for brown dwarfs, but the scattering 
results for the planet simulations are very similar. 
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Fig. 4. — Evolution of a planet/brown dwarf binary system from 
model t21 resulting in tidal capture. Conventions are as in Fig. 
[3] The system becomes dynamically active shortly after the calcu- 
lation begins. A strong tidal encounter at a later time binds the 
planet -planet pair, resulting in a configuration that is dynamically 
stable. The planet-planet binary's circularized semi-major axis is 
2.3 Rq. 



start of the calculation. Over time the planets exchange 
position relative to the star. The planets suffer a weak 
dissipative tidal encounter at late time but eventually 
collide and merge. 

Fig. [4] shows a typical "direct" tidal capture simula- 
tion. Tne system becomes dynamically active shortly 
after the calculation begins. A strong tidal encounter 
at a later time binds the planet-planet pair, resulting in 
a configuration that is dynamically stable. The planet- 
planet binary's circularized semi-major axis is 2.3 Rq. 
Fig. [5] on the other hand, shows a typical "gradual" tidal 
capture simulation. A weak tidal encounter (the smaller 



of the two open circles) nearly results in a dynamically 
stable planet-planet binary. The planets later suffer a 
stronger tidal encounter (the larger of the two open cir- 
cles) that results in a stable configuration. In this case, 
the planet-planet binary's circularized semi-major axis is 
6.0 i?0. These "gradual" tidal captures are responsible 
for the long tails in the final semi-major axis distributions 
in Figs. 1 and 2. 

2.6. Summary of Scattering Results 

For reference, we give the outcome statistics for each 
model in Table [2j We also show in Fig. [6] a summary 
of the fraction of dynamically active systems in which 
one of the planets is ejected as a function of the different 
assumed distances from the parent star, and in Fig. [Tithe 
fraction of systems leading to a successful tidal capture. 

2.7. Simulated Population of Binary Planets 

Based on the results of the numerous scattering events 
that we have computed for a wide range of masses and 
distances from the parent star, we have devised a sim- 
ple fitting formula that approximates the distribution 
for forming a tidally-captured binary pair of planets or 
brown dwarfs with a (final) circularized orbital separa- 
tion, flcirc- Expressed in Monte Carlo form, we have: 



:W 



R 



'i, bloat 



o.ix 
1 -X 



R, 



•0 I 



(3) 



where i?i, bloat is the thermally bloated radius of the larger 
of the planets/brown dwarfs at the time of tidal capture, 
i?j is the radius of Jupiter, and AT is a uniformly dis- 
tributed random number between and 0.99 (the upper 
limit in X introduces an upper cutoff of Ocirc — lOi?0, 
above which our scattering experiments provide insuffi- 
cient statistics). While the overall probabilities for tidal 
capture do vary systematically with distance from the 
parent star (see Table [2]) and the masses of the scatter- 
ing objects, the basic functional form of the distribution 
of Ocirc, and its dependence on the radii of the tidal cap- 
turing planets, seem relatively independent of distance 
and masses. 

Once we have in place analytic expressions for the ther- 
mal bloating as a function of age, and the distribution 
of tidal-capture circularized orbital radii as a function 
of the thermal bloating, we can use these to generate a 
synthetic population of binary planets and brown dwarfs. 
In this simulation, we first choose a random location, d, 
for the planet/brown dwarf binary between 0.01 and 10 
AU, uniformly distributed in log d. This is completely 
arbitrary since we do not know a priori the initial dis- 
tribution of where planets and brown dwarfs form in 
a protoplanetary disk. Next, we choose a random age 
for the planet/brown dwarf pair between 1 and 15 Myr, 
uniformly distributed in logt. The age sets the size of 
the thermally bloated radius according to equation (pi). 
The final (late-age) circularized orbital separation, Ccircj 
of the binary is chosen randomly from the distribution 
given by equation p]). Finally, we choose the eccentric- 
ity of the outer orbit (i.e., the orbit of the CM of the 
planet /brown-dwarf binary around the parent star) from 
the following probability distribution derived from our 
scattering studies: 



p{e) ex exp(— e/0.Q5) . 



(4) 



Podsiadlowski, et al. 




Fig. 5. — Evolution of a planet/brown dwarf binary system from 
model t21 resulting in a somewhat gradual tidal capture. Con- 
ventions are as in Fig. [s] A weak tidal encounter (the smaller of 
the two open circles) nearly results in a dynamically stable planet- 
planet binary. The planets later suffer a stronger tidal encounter 
(the larger of the two open circles) that results in a stable config- 
uration. Note that the radius of the open circles is proportional 
to the fractional change in relative planet-planet velocity result- 
ing from the encounter. In this case, the planet-planet binary's 
circularized semi-major axis is 6.0 Rq. 



We require for each binary pair, chosen according to 
the above prescriptions, that its orbital separation Ocirc 
be smaUer than 40% of the radius of the HiU sphere for 
that particular binary (this depends on the masses and 
the eccentricity of t he outer binary) in or der to ensure 



long-term stability ( Domingos et al. 2006 ) ^ Once the 
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Fig. 6. — The fraction of dynamically active planet systems in 
which one of the planets is ejected as a function of the distance 
from the parent star. Triangles are for "planets" in the brown- 
dwarf mass range, while circles represent gas giants of approxi- 
mately Jupiter's mass. Symbol size is proportional to the log of 
the planet/brown dwarf radius. 



ft 
O 





1 ' 




- 


o 

o 


6 
o 
o 


- A 


A 

A 


4 - 


A 



A 


o 
o 


A 


A 






- 


o 


- 


o 




- 


1 1 1 


1 






10-' 



Distance [AU] 



^'^ Their stability criteria apply to the case of massless satellites 
of planets. We have performed some numerical tests for equal- 
mass planet-planet binaries in circular orbits and obtained similar 



Fig. 7. — The fraction of dynamically active planet systems in 
which a successful tidal capture occurs leading to the formation of 
a stable binary planet. Triangles are for "planets" in the brown- 
dwarf mass range, while circles represent gas giants of approxi- 
mately Jupiter's mass. Symbol size is proportional to the log of 
the planet/brown dwarf radius. 



parameters of a binary have been fully chosen, we assign 
a probability of detection via transits simply as oc i?* /d, 



stability boundaries. 
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Fig. 8. — Simulated populations of Jovian binary planets (bot- 
tom panel), super- Jupiter binary planets (middle panel) and bi- 
nary brown dwarfs (top panel) in orbit about the parent star 
(in this case, IA/0). The masses of each binary pair are 
{70, 30}, {20, 10}, {1, 4} Jupiter masses (top to bottom panels). 
The color shading indicates the relative probability of discover- 
ing such a system, assuming that the binaries are born uniformly 
distributed in log d. White to purple colors indicate probability 
ratios of ~ 50 : 1. See text for specific algorithms used here. The 
radius of the Hill sphere is indicated by the diagonal green line. 



where d is the distance of the binary CM from the parent 
star of radius i?* . Obviously, this probability could be en- 
hanced if the projected separation of the planet/brown- 
dwarf binary, as seen from the Earth, is larger than the 
parent star. However, such calculations are beyond the 
scope of the paper. 



The results of our simulations are shown in Fig. [8] for 
three different mass pairs of planets/brown dwarfs. The 
top, middle, and bottom panels are for the mass pairs: 
{70, 30}, {20, 10}, {1, |} Jupiter masses, respectively. In 
each panel the color shading is proportional to the rel- 
ative probability of finding a transiting binary pair at 
planet distance, d, and planet separation, acirc- The ra- 
dius of the Hill sphere is shown as the diagonal green 
line. 
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Fig. 9. — The ^ 80 known transiting exoplanet systems (taken 
from "The Extrasolar Encyclopedia" http://exoplanet.eu/catalog^ 
all ) in the planet mass— semimajor axis plane, i'he blue line divides 
the allowed from excluded region of this space if the binary planet 
must fit within 40% of its Hill sphere and the binary components 
arc separated by at least 3 Rj. The dashed line is slightly more 
conservative (see text). 



In this work we have shown that if planets are dy- 
namically active early in their history, and that our pre- 
scription for tidal interactions is plausibly accurate, then 
binary planets should form with comparable frequency 
with which others are ejected. Therefore, since dynami- 
cally active planets are the currently favored me chanism 



for the formation of eccentric exoplanet orbits ( Ford h 
IRasio| [2008{ |Juric fc Tremaine| |2008| [Chatterjee eta 



2008), we expect that there exist gas-giant binary plan 
ets among the '-^ 420 that are currently known. In this 
section we discuss the signatures that such hypotheti- 
cal binary planets would exhibit in observations of the 
currently known sample of exoplancts, and prospects for 
detecting them in future studies. 

Since most of the known exoplanets have been discov- 
ered via radial velocity (hereafter, "RV") measurements, 
the first question to answer is how large would be the per- 
turbations to the observed Doppler signature. Treating 
the binary planet as an orbital perturbation, we derived 
an expression for the maximum deviations from a con- 
ventional RV curve of the central star due to the influence 
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TABLE 1 

Parameters for planet evolution simulations. 



model name 


mi/lO"^M0 


m2/10'^ Mq 


Ri/R@ 


R2/Rq 


ai,init/AU 


lidal dissipalion 


t2 


1 


1 


0.1 


0.1 


5 


off 


t3 


3 


3 


0.1 


0.1 


5 


off 


t4 


4 


2 


0.1 


0.1 


5 


off 


t5 


4 


2 


0.1 


0.1 


5 


on 


t6 


50 


25 


0.1 


0.1 


5 


on 


t7 


4 


2 


0.1 


0.1 


0.2 


on 


t8 


50 


25 


0.1 


0.1 


0.2 


on 


t9 


1 


0.5 


0.1 


0.1 


5 


on 


tio 


70 


30 


0.1 


0.1 


5 


on 


til 


1 


0.5 


0.1 


0.1 


1 


on 


tl2 


70 


30 


0.1 


0.1 


1 


on 


tl3 


1 


0.5 


0.1 


0.1 


0.2 


on 


tl4 


70 


30 


0.1 


0.1 


0.2 


on 


tl5 


70 


30 


0.2 


0.2 


0.2 


on 


tl6 


1 


0.5 


0.2 


0.2 


5 


on 


tl7 


70 


30 


0.2 


0.2 


5 


on 


tl8 


1 


0.5 


0.2 


0.2 


1 


on 


tl9 


70 


30 


0.2 


0.2 


1 


on 


120 


1 


0.5 


0.2 


0.2 


0.2 


on 


121 


70 


30 


0.6 


0.4 


0.2 


on 


122 


1 


0.5 


0.4 


0.4 


5 


on 


123 


70 


30 


0.4 


0.4 


5 


on 


124 


1 


0.5 


0.4 


0.4 


1 


on 


125 


70 


30 


0.4 


0.4 


1 


on 


126 


1 


0.5 


0.4 


0.4 


0.2 


on 


127 


70 


30 


0.4 


0.4 


0.2 


on 


128 


1 


0.5 


0.6 


0.6 


5 


on 


129 


70 


30 


0.6 


0.6 


5 


on 


130 


1 


0.5 


0.6 


0.6 


1 


on 


131 


70 


30 


0.6 


0.6 


1 


on 


132 


1 


0.5 


0.6 


0.6 


0.2 


on 


133 


70 


30 


0.6 


0.6 


0.2 


on 


134 


1 


0.5 


0.8 


0.8 


5 


on 


135 


70 


30 


0.8 


0.8 


5 


on 


136 


1 


0.5 


0.8 


0.8 


1 


on 


137 


70 


30 


0.8 


0.8 


1 


on 


138 


1 


0.5 


0.8 


0.8 


0.2 


on 


139 


70 


30 


0.8 


0.8 


0.2 


on 


140 


1 


0.5 


0.2 


0.2 


0.04 


on 



Note. — The quantilies nii and Ri are Ihe planel masses and radii, and aiinit is the initial 
semi-major axis ol planet 1. At least 5000 simulations were run lor each model. 



of a binary planet as 

9 
32 



For illustrative values oi m — Mj and Vb = 30 km s 



the numerical value of AVI' 



is ~ 0.3 cms , far too 



a 

Rq 



no 

UJ 



2m 
If 



^0, 



(5) 



where a is the orbital separation of the binary planets, 
Rq and JIq are the mean orbital radius of the outer binary 
and its mean angular velocity, respectively, Vq — QoRq, 
UJ is the synodic orbital period of the binary planet, m is 
the mass of an assumed equal-mass member of the binary 
planet, and M is the mass of the parent star. The mean 
orbits of the inner and outer binaries have been taken as 
circular and coplanar for simplicity. Since the ratio VLq/uj 
can be expressed in terms of the masses and orbital radii 
involved, we can write 



AK 




(6) 



For a binary planet with an orbital separation of a = 
0.4i?Hiii ^ 0.4(2to/3M)^/3^ ie., the largest stable sepa- 
ration, we have 



AK" 



9 
32 



3757 



7/6 



2m 

~M 



5/3 



K). 



(7) 



small to be detected in current observations as well as 
those that are currently planned for the near future. On 
the other hand, for a binary brown dwarf of, e.g., 50 Mj, 
the expected perturbations to the RV curve would be of 
order 2 m s^^, which would be detectable. Thus, it seems 
fair to say that the existence of gas-giant planets would 
not have been noticed in exoplanet RV curves. 

We next consider whether binary planets would have 
been detected among any of the ^ 80 currently known 
transiting exoplanets. For these systems, the transit light 
curves would be manifestly anomalous and the presence 
of two planets would be quite obvious. In Fig.|9]we show 
a plot of 80 currently known transiting exoplanets in the 
planet mass-semimajor axis plane. Note that, as ex- 
pected, most of these systems are close to the parent star 
(i.e., < 0.1 AU), thereby enhancing their transit proba- 
bility. The solid blue line indicates where in this plane a 
binary planet can fit within 40% of its Hill sphere and still 
allow the binary to be separated by at least 3 times the 
radius of Jupiter, thereby avoiding Roche lobe overflow 
from one planet to another. In this conservative set of 
restrictions, we see that only ^ 1/4 of the systems could 
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TABLE 2 
Outcome numbers for planet evolution simulations. 



model name 


total 


ejection 


collision 


tidal capture 


two planets 


stargrazer 


error^ 


t2 


5000 


1450 


2359 





1124 


6 


61 


t3 


5000 


2386 


1777 





830 


5 


2 


t4 


5000 


2436 


1628 





805 


123 


8 


t5 


5000 


1969 


1377 


969 


642 


33 


10 


t6 


5000 


3934 


373 


272 


395 


25 


1 


t7 


5000 


335 


3804 


286 


550 


25 





t8 


5000 


2760 


1426 


376 


303 


127 


8 


t9 


5000 


1234 


2112 


1337 


309 


1 


7 


tio 


5000 


4275 


301 


212 


165 


47 





til 


5000 


401 


3348 


1121 


127 





3 


tl2 


5000 


3905 


612 


313 


106 


64 





tl3 


5000 


119 


4653 


155 


71 


1 


1 


tl4 


5000 


3461 


1030 


320 


41 


148 





tl5 


5000 


3213 


1334 


302 


31 


119 


1 


tl6 


5000 


834 


2685 


1262 


214 


1 


4 


tl7 


5000 


4112 


428 


279 


144 


37 





tl8 


5000 


194 


4042 


666 


97 





1 


tl9 


5000 


3734 


771 


360 


78 


57 





t20 


5000 


82 


4809 


40 


67 


2 





t21 


100000 


58772 


36438 


2389 


532 


1856 


13 


t22 


5000 


495 


3204 


1166 


133 





2 


t23 


5000 


3957 


571 


317 


118 


37 





t24 


5000 


128 


4556 


238 


76 





2 


t25 


5000 


3517 


1001 


365 


55 


62 





t26 


5000 


43 


4892 


4 


61 








t27 


5000 


2959 


1743 


177 


27 


94 





t28 


5000 


369 


3548 


961 


120 





2 


t29 


5000 


3885 


657 


312 


104 


41 


1 


tso 


5000 


95 


4715 


122 


67 





1 


tsi 


5000 


3404 


1132 


357 


47 


60 





t32 


5000 


43 


4897 


3 


57 








t33 


5000 


2828 


1928 


134 


27 


82 


1 


t34 


5000 


265 


3851 


788 


94 


1 


1 


t35 


5000 


3801 


746 


325 


97 


30 


1 


t36 


5000 


97 


4791 


49 


62 


1 





t37 


5000 


3277 


1256 


358 


43 


63 


3 


t38 


5000 


46 


4898 


1 


54 


1 





t39 


5000 


2727 


2097 


67 


27 


82 





t40 


5000 


37 


4902 





60 


1 






Note. — Wc label the outcomes as inpord Sz Rasio (2008). "Two planets" refers to a system 
that has not achieved another outcome (collision, ejection, etc.) within the maximum integration 
time of 5 X 10® code units (^ 8 X 10^ initial orbital periods of planet 1). 

" For reference, we give the outcome statistics for each model. The final column includes both 
errors and uncounted outcomes. A typical error is the system becoming dynamically stable (and 
classified as such) due to integrator energy drift. A typical uncounted outcome is one planet being 
ejected while the other collides with the host star. 



possibly harbor a binary planet. If we make these con- 
straints only somewhat more 'comfortable' by requiring 
that the planets be separated by at least 5 Rj and that 
this separation be less than 1/5 of the radius of the Hill 
sphere, the separatrix is shown with a dashed blue line. 
In that case, it seems quite plausible that no more than 
three of the currently known transiting systems could 
contain a binary planet. At the moment, the numeri- 
cal simulations presented in this work do not place tight 
constraints on the ratio of highly eccentric exoplanet or- 
bits and those with binary planets, nor can they predict 
an absolute fraction of exoplanets that should be binary. 
Nonetheless, it is by no means obvious that any binary 
planets should have yet been detected by either RV or 
transit light curve studies. 

It is also possible that binary planets could be detected 
with gravitational lensing. The Einstein radius, Re, of a 
single planet of mass, m, at a distance D, and a generic 
source at ^ 2D, projected back to the lens plane has a 



physical size 



-Re 



2GMD 



0.06 AU 



mD] 



kpc 



M, 



(8) 



Thus, if the planet is actually a binary, it would have to 
have an orbital separation comparable to, or greater than 
Re, in order for its binary nature to be readily revealed 
in the microlensing light curves. For comparison, the 
maximum separation between binary components of a 
Jupiter mass at an AU from a 1 Mq star is ~ 0.03 AU. 

Probably the best hope for detecting binary planets is 
to obtain a larger sample of transiting exoplanets, yield- 
ing a substantial number (e.g., > a dozen) at distances 
from the parent star of > 0.4 AU. The Kepler mission 
should ultimately provide such a sample. At these larger 
distances for transiting exoplanets, there is more phase 
space to fit stable binaries within their Hill sphere and to 
allow larger separations between the binary components. 
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Fig. 10. — Left Panel: Comparison with Fig. 2 of [Ford &: Rasio | ( |2008|l , the cumulative e ccentricity distribution of the remaining planet 
ifter a planet is ejected for the case mi /Mi, = m^lMi, = 10~'^. As in|fbrd fc Rasiol 1120081), we set the stellar mass to Mi, = 1 Mq and the 



after a planet is ejected for the case mi/Mi, = m^lMi, 

planet masses to mi = 10~^ ^Q- Right Panel: Comparison with Fig. 3 of |Ford fc Rasio ( 2008 ) , the cumula tive eccentricity distribution of 

the remaining planet after a planet is ejected for /3 ; 

and Mi, = IMq. 



mi/m2 = 1/2 and = 1/3. As in |Ford fc Rasio| ( |2008t , we set mi +7712 = 6 X 10"^M* 



If gas-giant binary planets are ultimately discovered, 
they would represent a strong corroboration of the dy- 
namically active scenario, since gas giant binaries formed 
via tidal capture are a natural outcome of dynamically 
active systems. On the other hand, if no binary plan- 
ets are detected, even with a larger sample of transiting 
exoplanets, it may simply be that our treatment of the 
tidal capture process is too optimistic. Parameters of 



the tidal capture process could be constrained by com- 
paring A^-body simulations that allow for tidal capture 
with observed planet eccentricity distributions, with the 
restriction that no observable planet-planet binaries be 
formed. 



JMF acknowledges support from Chandra Postdoc- 
toral Fellowship Award PF7-80047. 



APPENDIX 

USE OF THE FEWBODY DYNAMICAL CODE 

In our nu merical simu lations we use the Fewbody integrator, which is designed for strong small- A'^-body gravitational 
encounters ( Fregeau et al._^2004j . To test Fewbody^ s sui tability for dynamically active planetary systems, we compare 
our scattering calculations with the work of Ford & Rasio (2008), who studied planet-planet scattering with the aim of 
explaining the high eccentricities of some observed extrasolar planetary systems. As with mo st studies of planet-planet 
scattering they used a mixed variable symplectic method modified to treat close encounters ( Wisdom fc Holman|1991 
Chambers||1999 1 . When this method is applied to a two-planet system, a close encounter between the planets results 



in all orbital motion being integrated with a standard, non-symplectic integrator (e.g., Bulirsch-Stoer). Since the two- 
planet systems we study are dynamically active and hence quickly result in close approaches, their evolution should 
be fai thfu lly treated (in a statistical sense) with the adaptive, but non-symplect ic, integration algorit hm in Fewbody. 



Fig. [To] (left panel) shows a comparison of our numerical method with Fig. 2 of [Ford fc Rasio ( 2008 1 , who integrated 
the evolution of dynamically active two-planet systems. Specifically, the cumulative eccentricity distribution of the 
remaining planet after a planet is ejected is shown for the case mi /Mi, — mijMi, = 10"'^, where nii is the mass of planet 
i and Mi, is the mass of the host star. By convention, planet 1 initially orbits with semi-major axis ai^injt < a24nit- Our 
initial conditions ( including planet ecce ntricities, relative inclinations, and randomization of various orbital elements) 
are the same as in Ford & Rasio (2008). Most notably for this case, we similarly set the stellar mass to Mi, — \Mq 
and the planet masses to m^ = 10"'^ M©. As is clear from the figure, the eccentricity distributions are fairly similar - 
in fact, the median eccentricity is the same in both distributions. The differences are most ap parent at low and high 
eccentricities. Part of the discrepancy may arise from our differing definitions of ejection. In Ford &: Rasio (12008), 
a planet straying beyond 2000ai.init is declared to have been ejected, while in Fewbody a planet is not considered 
ejected until it is formally unbound from and receding from the remaining masses, and is sufficiently weakly tidally 
coupled to the remaining masses that it could never become bound again. Another source for t he discrepancy may b e 
the accuracy of th e mixed variable symplec tic method for close star-planet approaches used in Ford & Rasio (2008). 
As pointed out in iRauch & Holman (1999), the method introduces artificial chaos for close star-planet approaches. 
Finally, a remaining source tor the discrepancy may be the calculation stopping time. We used a stopping time of 
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Ford & Rasio (2008) used a 



5 X 10^ code units (^ 8 x 10^ initial orbital periods of planet 1) for all calculations 
stopping time between 5 x 10® and 2 x 10^ that is an unspecified func tion of the planet ma sses 

Fig. 10 (right panel) shows a comparison of our results with Fig. 3 of Ford & Rasio (20081. Shown is the cumulative 
eccentricity distribu tion of the remaining planet after a planet is ejected for jj = mi/m2 = 1/2 and /3 = 1/3. As in 
Rjrdfc Rasio1(|2008|), we set M^ 



6 X 10 '^M^ to speed up the evolution. (As shown in [Ford fc 
Rasio (2008), the final eccentricity distribution is not sensitive to the ratio (mi +m2)/M* for mi/ Mi, ~ 10""^, but the 



1 Mq and mi + ra2 



time to ejection decreases as (mi + m2)/M* increases.) While the agr e ement between the eccentricity distributions 
is not perfect, our results agree fairly well with those of Ford & Rasio ( 2008| ), and reproduc e the dependence of t he 
eccentricity distribution on /3. Additionally, eccentricities above 0.8 are very rare, as found in |Ford fc Rasio (2008). 
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